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Abstract 

We analyzed the photoluminescence intermittency generated by a single paramagnetic spin lo- 
calized in an individual semiconductor quantum dot. The statistics of the photons emitted by the 
quantum dot reflect the quantum fluctuations of the localized spin interacting with the injected car- 
riers. Photon correlation measurements which are reported here reveal unique signatures of these 
fluctuations. A phenomenological model is proposed to quantitatively describe these observations, 
allowing a measurement of the spin dynamics of an individual magnetic atom at zero magnetic 
fleld. These results demonstrate the existence of an efficient spin relaxation channel arising from a 
spin-exchange with individual carriers surrounding the quantum dot. A theoretical description of 
a spin-flip mechanism involving spin exchange with surrounding carriers gives relaxation times in 
good agreement with the measured dynamics. 

PACS numbers: 78.67.Hc, 78.55.Et, 75.75.+a 



1 



I. INTRODUCTION 



The decrease of the structure size in semiconductor electronic devices and magnetic infor- 
mation storage devices has dramatically reduced the number of atoms necessary to process 
and store bit of information. Information storage on a single magnetic atom would be an 
ultimate limit. The performance of such memory elements will be governed by the quan- 
tum fluctuations of the localized spins^. Diluted magnetic semiconductors (DMS) systems 
combining high quality semiconductor structures and the magnetic properties of Mn impu- 
rities are good candidates for these ultimate spintronic devices- . It has been shown that in 
a DMS with low magnetic atom concentration, the spin dynamics under magnetic field is 
ultimately controlled by the spin-lattice coupling^i^. An extrapolation of the spin dynamics 
measurements in bulk DMS suggests that a long spin relaxation time in the millisecond 
range could be expected for an isolated Mn spin^. However, despite the recent development 
of different successful techniques to address a single spin^i^, such dynamics has never been 
directly observed. 

Growth and optical addressing of DMS quantum dots (QDs) containing a single magnetic 
atom were achieved recently^i^. When a Mn dopant atom is included in a II- VI QD, the 
spin of the optically created electron-hole pair interacts with the five d electrons of the Mn 
(total spin M=5/2). This leads to a splitting of the once simple photoluminescence (PL) 
spectrum of an individual QD into six (2M+1) components. This splitting results from the 
spin structure of the confined holes which are quantized along the QDs' growth axis with 
their spin component taking only the values J2=±3/2. The hole-Mn exchange interaction 
reduces to an Ising term J^.M^ and shifts the emission energy of the QD, depending on the 
relative projection of the Mn and hole spins. As the spin state of the Mn atom fluctuates 
during the optical measurements, the six lines are observed simultaneously in time average 
PL spectrum. The intensities of the lines reflect the probability for the Mn to be in one of 
its six spin components and the emitted photon is a probe of the spin state of the Mn when 
the exciton recombines. 

In this article, we show how we can use the statistics of the photons emitted by a single 
Mn doped QD to probe the spin dynamics of a magnetic atom interacting with both confined 
and free carriers. We performed correlations of the PL intensity emitted by individual lines 
of an isolated Mn-doped QD. In these start-stop experiments^, the detection of the first 
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FIG. 1: (Color online) (a) Experimental time resolved PL spectra recorded on three different 
energy lines (labeled 1, 2 and 3) of an X-Mn complex. The inset shows the corresponding time 
integrated PL spectrum, (b) PL decay time calculated using the parameters T=5K, r;,=280ps, 
T(i=8ns, Tx-Afn=25ns. (c) Energy levels involved in the rate equation model described in the text 
displayed as a function of their total angular momentum S^. 

photon indicates by its energy and polarisation that the Mn spin has a given orientation. 
The detection probability of a second photon with the same energy and polarisation is 
then proportional to the probability of conserving this spin state. The time evolution of 
this intensity correlation signal is a probe of the spin dynamics of the Mn atom. The auto- 
correlation signal displays a bunching effect revealing a PL intermittency. This intermittency 
results from fluctuations of the Mn spin. Correlation of the intensity emitted by opposite 
spin states of the exciton-Mn complex (namely,cross-correlation) presents an antibunching 
at short delays. The characteristic time of this antibunching corresponds to the spin transfer 
time between the two degenerated spin states. A thermalization on the exciton-Mn complex 
is directly evidence by the energy and temperature dependences of the correlation curves. 
The measured single Mn spin relaxation times are in the range of 20 ns and are strongly 
influenced by the injection of carriers in the vicinity of the QDs. Mn spin-flips induced by 
the injection of carriers in and around the QD are theoretically described. These scattering 
processes give relaxation times in good agreement with the measured dynamics. 
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II. SAMPLES AND EXPERIMENTS 



Single Mn atoms are introduced in CdTe/ZnTe QDs during their growth by molecular 
beam epitaxy adjusting the density of Mn atoms to be roughly equal to the density of QDs^. 
The statistics of the photons emitted by individual Mn-doped QDs was analyzed using a 
combination of a low-temperature (5K) scanning optical microscope and a Hanbury Brown 
and Twiss (HBT) setup for photon-correlation measurements^S,. High refractive index hemi- 
spherical solid immersion lens were used to enhance the collection of the single dot emission. 
The PL was quasi-resonantly excited with a tunable CW dye laser or by picosecond pulses 
from a doubled optical parametric oscillator and collected through aluminum shadow masks. 
The circularly polarized collected light was spectrally dispersed by a 1 m monochromator 
before being detected in the HBT setup or by a single fast avalanche photodiode (time 
resolution 40 ps) for time resolved measurements. The time resolution of the HBT setup 
was about 500 ps. Under our experimental conditions with counts rates of a few kHz the 
measured photon pair distribution yields after normalization the intensity autocorrelation 
function g*^^)(r). 

III. EXCITON-MN SPIN FLIPS. 

A signature of the Mn spin dynamics can be observed first in the PL decay of the X-Mn 
complex. Fig.l presents the PL decay of three different transitions of the X-Mn complex. 
These transitions present a biexponential decay. The fast component corresponds to the 
radiative lifetime of the exciton, as already measured in non magnetic QDs. The long 
component is associated with the existence of the dark excitons^^. Direct recombination 
of the dark exciton can be observed in some Mn-doped QDs because of a slight admixture 
of the bright states with the dark ones induced by a valence band mixing. However, the 
dark excitons mainly contribute to the signal by undergoing a spin flip to become bright 
excitons which decays radiatively. The PL decay is then determined by both radiative 
decay and excitons spin-flips. The exciton decay, and particularly the amplitude of the slow 
component, depends strongly on the energy level observed. For the high energy lines, the 
slow component makes a negligible contribution to the time integrated signal. Conversely, 
for the low energy lines, the secondary component makes a significant contribution while the 
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FIG. 2: (a) Auto-correlation function of the intensity collected in a+ polarisation on the low energy 
line of the X-Mn complex (solid line) and on the overall PL spectra (dotted line), (b) Circularly 
polarized cross-correlation function (solid line) and auto-correlation (dotted line) on the same line 
as in (a) but for a larger excitation intensity. The inset shows the PL spectrum of the corresponding 
QD. (c) Experimental auto-correlation function and theoretical function calculated with the rate 
equation model described in the text with parameters T=5K, rfe=280ps, Td=8ns, Tx-Afn=25ns, 
TMn=50ns and g=0.05xl0^^ps^^. 

primary lifetime remains constant. In this regime, the secondary lifetime can be associated 
either with the dark exciton lifetime or the exciton spin flip time. 

To extract these two parameters from the PL decay curves, we compare the experimental 
data with a rate equation model describing the time evolution of the population of the 24 
X-Mn spin levels (Fig. 1(c)) after the injection of a single exciton^^. Different spin-flips times 
are expected depending on wether the transitions occur with or without conservation of the 
energy or of the total spin. However, we consider in first approximation that the spin-flips 
among the X-Mn states can be described by a single characteristic time tx-mu- We consider 
that at finite temperature, the intraband relaxation rates T^^^i between the different spin 
states of the exciton- Mn complex depend on the energy separation E^^fi = Ey —E^. Here we 
use T^^Y = 1/Tx-Mn if E^y < and T^^y = l/Tx-Mne~^'''''^''^^ if E^y > 0^^- This describes 
a partial thermalization among the 24 levels of the X-Mn system during the lifetime of the 
exciton (bright or dark). In this model, we also neglect the influence of the valence band 
mixing on the oscillator strengtb^^ and consider that all the excitonic bright (dark) states 
have the same lifetime Th (r^). Because of an efficient hole spin-flip during the phonon 
assisted relaxation of the unpolarized injected carriers, we consider that the excitons with 
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spins ±1 and ±2 are excited with the same probabihty at t = 0. Only optical transitions 
conserving the Mn spin are considered. The calculation shows that changing mainly 
influences the characteristic time of the long decay component whereas changing tx-mu 
affects also the amplitude of this component. The PL decay curves can be reproduced well 
by this rate equation model using r;, = 280ps, = 8ns and tx-mu = 25?2S (Fig. 1(b)). The 
value of Td controls the decay time of the long component whereas tx-mu, larger than r^, 
reproduces the emission energy dependence of the amplitude of the long component very 
well. 

Within the X-Mn complex, the relaxation time between different spin states tx-mu can be 
affected by the spin orbit, carrier-phonon and exchange interactions— that affect the exciton. 
In QDs, the electron spin relaxation is longer than the radiative lifetime and is ultimately 
controlled by random fluctuations of the nuclear spins^^. The hole spin relaxation is mainly 
controlled by the interaction with phonons^i and can be faster explaining the observed 
partial thermalization of the X-Mn complex during the lifetime of the ground state exciton. 

IV. FLUCTUATIONS OF AN ISOLATED MN SPIN. 

To directly observe the time fluctuations of the Mn spin interacting with the injected 
carriers, we analyzed the statistics of the photons emitted by a Mn-doped QD. This statistics 
can be deduced from an intensity correlation function of the QD emission. Fig. 2(a) shows 
the intensity correlation function g^^-'('r) of the circularly polarized (o"+) low energy line of a 
Mn-doped QD compared with the correlation function obtained for the overall PL of the QD. 
The auto-correlation function obtained for all the photons emitted by the QD is characteristic 
of a single photon emitter with a dip at short delay. The width of this antibunching signal 
is given by the lifetime of the emitter and the generation rate of excitons and its depth is 
limited by the time resolution of the HBT setup. A similar experiment performed on one 
of the single line of the X-Mn complex still presents a reduced coincidence rate near r=0, 
but it is mainly characterized by a large bunching signal with a half width at half maximum 
(HWHM) of about 10ns. This bunching reflects an intermittency in the QD emission. This 
intermittency likely comes from the fluctuations of the Mn spin orientation. 

To confirm this result, cross-correlation measurements were performed between different 
spin states of the X-Mn complex. Cross-correlation of the cr+ and cr— photon emitted 
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FIG. 3: (Color online) (a) Power dependence of the autocorrelation function of the hight energy 
line of an X-Mn complex, (b) Calculated power dependence of the HWHM of the autocorrelation 
function of the hight energy line. The parameters used in the model are: r{,=280ps, Td,=8rLS, 
Tx-Mn=20ns, and rMn=40ns (plain line) or TMn=^IJ-^ (dotted line), (c) Experimental HWHM of 
the autocorrelation signals presented in (a). 

by the low energy line (fig.2(b)) shows an antibunching with g(^)(0)=0.2 and a HWHM of 
about 5ns. These two different behaviors, namely the bunching of the auto-correlation signal 
and the antibunching of the cross-correlation signal, demonstrate unambiguously that the 
statistic of the QD emission is completely governed by the Mn spin fluctuations. Whereas 
the auto-correlation probes the time dependence of the probability for the spin of the Mn 
to be conserved (M^ —-\-b/2 at t=0 in Fig.2(a)), the cross-correlation presented in Fig. 2(b) 
is a probe of the spin transfer between -1-5/2 and -5/2. 

As observed in the time resolved PL measurements, fluctuations of the Mn spin occur 
during the lifetime of an exciton. However, they can also take place when the QD is empty. 
As the spin relaxation rate of the Mn is expected to be influenced by the presence of carriers 
in the QD, we have to consider two relaxation times, tmu for an empty dot and tx-mu 
for a dot occupied by an exciton. Their relative contributions to the observed effective 
relaxation time will depend on the generation rate of excitons. The rate equation model 
described previously can be extended to extract an order of magnitude of the parameters 
TMn and Tx-Mn from the correlation experiments. Six biexciton states are added to the 24 
{X — Mn)+Q {Mn) level system (see Fig. 1(c)). A continuous generation rate g is considered 
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to populate the exciton and biexciton states. The initial state of the system is fixed on one 
of the six Mn ground states and one monitors the time evolution of the population of the 
corresponding bright X-Mn state. When normalized to one at long time, this time evolution 
accounts for the correlation function of the transition associated with the considered X-Mn 
level. 

The time evolution of the correlation function calculated using this model are presented in 
Fig. 2(c) and compared with the experimental data. At low generation rate, when the average 
time between two injected excitons is longer than any spin relaxation rate, tmh and tx-mu 
have distinguishable effects on the calculated correlation curves. In average, the relaxation 
of the Mn alone (controlled by tmu) is only perturbed by the injection of the exciton used to 
probe the Mn spin projection. During the lifetime of this exciton, the system relaxes with 
the relaxation rate Tx-Mn- This produces a relaxation of the Mn spin proportional to the 
ratio of tx-mu and the exciton lifetime. A reduction of Tx-Mn then reduces the amplitude 
of the bunching curve expected for a Mn alone (because of the six available spin states, 
the maximum amplitude of the bunching should be six) without significantly changing its 
width controlled by tmu- With the generation rate used in the measurements of Fig. 2 (a) 
(a generation rate of about g=0.05xl0^'^ps^^ can be deduced from the ratio of the exciton 
and biexciton amplitude), tx-mu mainly affects the height of the bunching signal whereas 
TMn preferentially controls its width. Then, with given values of g, and t^, it is possible 
to extract a parameter pair {rMm'Tx-Mn) that reproduces the bunching and anti-bunching 
curves. The bright and dark excitons lifetimes were estimated from the PL decay curves 
and the exciton generation rate can be estimated from the relative amplitudes of the exciton 
and biexciton emissionsi^. For the data presented in Fig.2(c), the best fit is obtained with 
'Tx-Mn = '25ns and TMn = 50ns. The X-Mn relaxation time obtained is consistent with 
the value deduced from the PL decay curves. The relaxation time of the Mn alone (empty 
dot) appears to be 3 orders of magnitude shorter than expected from the extrapolation of 
measurements in bulk dilute CdMnTe under magnetic field^. 

V. CARRIERS INDUCED MN SPIN FLUCTUATIONS 

The observed Mn spin dynamics is not simply an intrinsic property of the localized 
Mn atom but depends on the optical excitation conditions. The power dependence of the 
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FIG. 4: (color online) (a) Auto-correlation function on the low energy line of an X-Mn complex 
in a+ polarization for excitation intensities Pq and SPq- Theoretical curves are presented in red. 
A reduction of Tx-Mn{=^5ns) and TMn(=20ns) has to be included to describe the high excitation 
intensity autocorrelation curve, (b) Auto-correlation function on the low energy line of an X-Mn 
complex in a+ polarization for two different excitation conditions: resonant on an excited state 
(577. 5nm) and non-resonant (514nm). 

correlation signal of the high energy transition of a X-Mn complex is presented in Fig. 3 (a). 
Increasing the excitation power significantly reduces the width of the correlation signal. This 
reduction has two origins: first, when carriers are injected in the QD under quasi-resonant 
conditions (excitation on an excited state of the QD), increasing the carrier generation rate 
increases the probability of finding the QD occupied by an exciton. The spin relaxation 
time being shorter for an occupied dot than for an empty dot, the Mn spin fluctuates faster 
and the width of the auto-correlation curve decreases. This effect is illustrated by the power 
dependence of the HWHM of the calculated and experimental correlation curves presented in 
Fig. 3(b) and 3(c) respectively. At high generation rate, the width of the correlation signal is 
controlled by tx^mu whereas at low excitation the photon statistics is ultimately determine 
by the spin fluctuations of the Mn alone. The width of the calculated correlation curves 
saturates at low excitation. This maximum width is controlled by tmu- In the experiments, 
this saturation is not observed due to the limit in the accessible excitation power range. 

However, this process is not sufficient to explain the observed excitation power dependence 
of the correlation signal in all the investigated QDs. For instance, to reproduce the power 
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dependence presented in Fig.4(a), one has also to reduce the spin relaxation times tmu 
and Tx-Mn at high excitation intensity. In Fig.4(a) the best fit at high excitation power 
is obtained with tmu = 25 and Tx-Mn = 15 whereas at low excitation tmu = 50 and 
'Tx-Mn = 25. This reduction of the relaxation time likely comes from the spin-spin coupling 
with carriers injected in the surroundings of the QD thought the background absorption 
observed in PLE spectra of this individual QI^. 

The infiuence of the free carriers on the spin relaxation rate is shown by the correla- 
tion signals obtained on the same X-Mn transition for two different excitation wavelengths 
(Fig. 4(b)): resonant on an excited state (577,5nm) and non-resonant in the ZnTe barriers 
(514nm). These two signals are recorded with the same photon count rate, suggesting a 
similar occupation rate of the QD. The characteristic bunching signal observed under quasi- 
resonant excitation completely disappears when the excitation energy is tuned above the 
wetting layer absorption. As already observed in DMS quantum wells, this behavior re- 
flects the extreme sensitivity of the localized Mn spin to the spin-spin coupling with the free 
carriers or the carriers relaxing in the QD^. 

For an isolated Mn atom, the spin relaxation tmu comes only from the spin-lattice 
interaction?^ and a long spin relaxation time is expected. As we will discuss in the next 
section, the Mn spin dynamics can be modifled signiflcantly by the presence of free carri- 
ers which are strongly coupled with both the magnetic atom and the phonons. These free 
carriers serve as a bypass channel for the slow direct spin-lattice relaxation. 

VI. THERMALIZATION OF THE EXCITON-MN COMPLEX. 

The X-Mn complex is also signiflcantly coupled to the phonon bath. A partial thermaliza- 
tion of the X-Mn system appears directly in the amplitude of the correlation curves obtained 
on different energy levels of the X-Mn system (Fig. 5 (a)) as well as in cross-correlation mea- 
surements (Fig. 5(b) and (c)). A flnite temperature enhances the probability of the spin-flips 
involving an energy loss. This introduces a dissymmetry in the spin relaxation channels 
of the X-Mn complex. The consequence of this dissymmetry in the relaxation process is 
an energy dependence of the amplitude of the correlation signal. This is illustrated by the 
correlation curves obtained on the high (| — 1, —5/2)) and low energy states (| + 1, —5/2)) of 
a X-Mn complex (Fig. 5 (a)). For the high energy state, all the relaxation transitions within 
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FIG. 5: (a) Auto-correlation function of the emission intensity of the hight (upper trace) and 
low (lower trace) energy lines of a X-Mn complex recorded in the same circular polarization, (b) 
Cross-correlation function of the emission intensity of the high and low energy line recorded in a+ 
and cr— polarization respectively for two different excitation intensities. Detail of the experimental 
(c) and calculated (d) cross-correlation function. 

the X-Mn complex take place with an energy loss: The leakage probability is then maximum 
and the probability for this state to be re-populated by spin-flips from low energy states is 
very weak. A large bunching signal is then observed (Fig.5(a)). On the opposite, the low 
energy level can be populated by a transfer from the high energy states, and some relaxation 
channels involving an absorption of energy are blocked at low temperature. The associated 
bunching signal is weaker (Fig. 5 (a)). 

This tlicrmalization process directly appears if a cross-correlation of the intensity emit- 
ted by a low and a high energy levels is performed. Fig. 5(c) shows the correlation of the 
photons emitted by | — 1,-5/2) (high energy line) and | -|- 1,-5/2) (low energy line). At 
low excitation intensity, this correlation signal presents a clear dissymmetry. This cross- 
correlation measurement probe the time dependence of the probability of finding an exciton 
(either | -|- 1) or | — 1)) coupled with the Mn spin in the state Mz—-5/2. At positive time 
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FIG. 6: Intensity auto-correlation function of the high energy hne | — 1, —5/2) of a X-Mn complex 
recorded in the same excitation conditions at T=5K (a) and T=20K (b). The theoretical curves 
(solid line) are obtained with the same set of parameters: Tx-Mn=Wns, rA,fn=15ns,rfe=280ps, 
rrf=8ns and g=0.1.10^^ps^^. 

delay, g^(T) gives the probability to find the system in the state | + 1, —5/2) knowing that 
at r = the Mn spin projection was Mz=-5/2 (detection of a photon from | — 1, —5/2)). 
The situation is reversed for the negative delay where a photon from | + 1, —5/2) acts as the 
trigger in the start-stop measurement and g^(r) give the probability for the system to be in 
the high energy state | — 1, —5/2). The dissymmetry in the cross-correlation curve reflects 
the difference in the spin relaxation channels available for the high (| — 1, —5/2)) and the 
low (I + 1, —5/2)) energy X-Mn states. 

The dissymmetry in the relaxation processes is influenced by the excitation intensity: as 
observed in the PL decay measurements presented in Fig.l, the low energy bright exciton 
states can be efficiently populated by spin-flips from the dark exciton states reducing the 
effective population loss of these states and consequently reducing the amplitude of the 
photon bunching. Increasing the excitation intensity decreases the effective lifetime of the 
dark excitons because of the formation of the biexciton™. This opens an efficient spin 
relaxation channel for the low energy bright X-Mn states: once a dark exciton has been 
created after a spin-flip, it is quickly destroyed by the injection of a second exciton with the 
formation of a biexciton. It can no longer flip back to the low energy bright state. This 
effect stop the refllling process and consequently increases the amplitude of the bunching 
signal. As observed in Fig.5(b), increasing the excitation intensity decreases the difference 
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in the amplitude of the corresponding bunching signal of the low and high energy lines. 

As shown in Fig. 6, these spin relaxation processes are also influenced by the lattice 
temperature. The effect is especially pronounced for the high energy states. Increasing the 
temperature allows a refilling of these levels by a transfer of population from the low energy 
states, decreasing the amplitude of the bunching signal. 

VII. MODEL FOR THE MN SPIN RELAXATION 

In this section we discuss two Mn spin relaxation mechanisms in the optical ground 
state. The short {tmu — 20ns) spin relaxation time inferred from our photon correlation 
measurements can not be accounted for by Mn-Mn spin interaction, which is considered the 
most efficient spin relaxation mechanism in II- VI semiconductor. Since this is a short range 
interaction—, the Mn spin relaxation time, Tim increases exponentially as the density of Mn 
goes down^. Thus, the spin relaxation of a single Mn atom inside a quantum dot with 20 
10^ atoms, if generated by Mn-Mn coupling, should be in the range of microseconds. Thus, 
some other spin relaxation mechanism is at play. Here we study the Mn spin relaxation 
due to exchange coupling to photo-carriers that occupy extended states in the wetting layer 
(WL). The Mn is assumed to lie in the wetting layer, at a location where both quantum 
dot and extended states are finite. In particular, the QD states are expected to peak in the 
plane that separates the wetting layer and the QD. Mn atoms located in that region are 
more strongly coupled to the QD exciton. 

We consider two independent mechanisms: (i)exchange coupling with carriers in extended 
WL states and (ii) exchange coupling with carriers that scatter from the extended WL 
states to confined QD states, exchange with the Mn and then return to the WL. The first 
mechanism has been considered before in the context of Mn spin relaxation interacting with 
quantum well carriers^^"^ and is identical to Korringa mechanism for nuclear spin relaxation 
due to spin-fiip with itinerant electrons in metals. The second mechanism involves single 
charge tunneling in and out of the dot. The QD confined states remains empty except for 
some intervals during which a single carrier tunnels in and out from the optically excited 
wetting layer. Once in the QD state, the carrier can give (or take away) one unit of spin to 
Mn. Both the Korringa and the charge fiuctuation mechanisms require several tunneling and 
exchange events, involving different carriers, to relax the Mn spin from an initial situation 
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where is known to take the value Mz = +5/2 to a final situation of thermal equilibrium, 
where the 6 Mn spin orientations are equally likely. 

The proposed charge fiuctuation induced spin relaxation mechanism is in line with our 
previous works^i^i22;^i22 whose central claim is the very strong dependence of the spin prop- 
erties of Mn doped quantum dots on the addition of a single carrier. It is also motivated by 
the observation of a weak contribution of the positive trion PL signal which implies that, 
under our experimental conditions, holes are captured by the QD. 

A. Formalism 

We use a Bloch-Redfield master equation approaches which tracks the dynamical evolu- 
tion of the density matrix of the quantum dot ground states, corresponding to the 6 Mn 
spin orientations, under the infiuence of the reservoir of carriers. The diagonal terms of the 
density matrix represent the probability of finding the Mn spin in a given spin projection. 
Their dynammics is given by the master equation 



with the initial condition that, at t = 0, the probability of finding the Mn spin with = 
+5/2 is p+5/2 = 1- The evolution of p+5/2(t) and p_5/2(t) are directly associated to the auto- 
correlation and cross-correlation measurements. The rates in the master equation depend 
on the Hamiltonian of the sytem, which is the sum of the QD Hamiltonian, the WL carrier 
Hamiltonian and the QD-WL coupling: H = Hq + Tires + V 

Tm^m' is the scattering rate between the eigenstate M and state M' of TIq induced by 
the coupling V. The rates are given by the statistical average over initial reservoir states 
and sum over final reservoir states of the Fermi Golden rule rate associated to V"^*^. They 
depend on the mechanism under consideration, either direct exchange or charge fiuctuation. 
Since our experimental data indicate that hole confinement is much weaker than electron 
confinement, we consider that the carriers that couple to the Mn when there is no exciton 
in the dot are holes. Our theory results can easily be adapted for the case of electrons. The 
number of states M in the master equation also depends on the mechanism. In the case 
of charge fiuctuation mechanism we need to keep track of the 6 neutral states and the 12 
states of the dot with one hole inside. In the case of the Korringa relaxation, only the 6 




(1) 
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states corresponding to the "empty dot" Mn spin are included in the master equation. 
The QD Hamiltonian is: 

710 = .] (r^m' + 1^ (t+MT + r-M+)) + D {yT^ (2) 

where J = l3\tpQD{ri)\'^nii is the Mn QD hole coupling, j3 ~ 60e^A^ is the Mn-hole exchange 
coupling constant of CdTe, \ipQD{r)\'^ is the QD hole envelope function, = 0, 1 counts 
the number of holes in the dot. The Mn spin operators are M^'± and r^'^ are the Pauli 
matrices acting on the hole space. We consider antiferromagnetic hole-Mn coupling, eh is 
the dimensionless parameter that accounts for the reduction of the spin-flip interaction due 
to spin-orbit coupling^i^^. 

Notice that the exchange coupling of the Mn to the quantum dot fermion is only relevant 
in the charge-fluctuation mechanism, for which Ud changes between and 1. For the Korringa 
mechanism Ud = and Mn is only exchanged coupled to carriers in extended WL states. 
The Z}(M^)^ term in Hqdo describes the strain induced anisotropy which has been observed 
in strained (Cd,Mn)Te layers^^. This term is negligible when the Mn interacts with a QD 
hole but is not when we consider the Korringa relaxation, for which no carrier occupies 
the QD state. The term Tires = J2k v ^k.i^cj^ yC\^,v describes the delocalized carriers in states 
momentum k and band index (including the spin) v. 

Each of the two relaxation mechanisms considered here has its own dot-reservoir coupling. 
In the Korringa mechanism we assume that there is some overlap between the WL extended 
states and the Mn spin. Igoniring the dependence on k of the spin matrix element e(k, k'),the 
coupling reads: 

V ^ E (M^+^i,Ck^^ + M(-)4,,Ck,,) (3) 

k,k' 

where A is the area of the WL, is the envelope part of the WL wave function 

evaluated at the Mn location, and L is the WL width, could take different values for the 
QD and the WL. 

In the charge-fluctuation mechanism the coupling between the reservoir and the dot is 
the sum of an operator that transfers one hole from the reservoir to the dot and its hermitian 
conjugate, which takes the hole out from the dot and transfers it to the reservoir: 

l^ = 5Zl^)(0|5Z^kCk,. + /i.c. (4) 

(J k 
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The tunneling operator conserves the energy and the spin of the carrier. This kind of couphng 
has been considered recently to account for peculiar PL lineshapes of self assembled quantum 
dots in contact with electronic reservoirs^^i^^. The mechanism would operate analogously if 
the transfer from the itinerant to the localized state is inelastic, but the timescale would be 
longer. 



B. Results 



1. Mn spin relaxation due to exchange with WL carriers 

We discuss first the relaxation of the Mn spin due to exchange with WL carriers, or 
Korringa mechanism. We assume a parabolic dispersion with effective mass m* for the holes 
which yields a stepwise density of states. We assume holes are in a thermal state with 
effective temperature ksT^ larger than the lattice and chemical potential /i. The Korringa 
relaxation rate is proportional to the square of the density of states at the Fermi energy and 
the square of the exchange coupling constant. After some work we obtain the Mn spin flip 
rates between M and M', eigenstates of Hq = -D(M^)^, : 

rt,^M' = (^)'/(y,^)|(M|MW|M')r (5) 

where we have approximated = 2r7/L,0<?7<lisa dimensionless parameter that 

accounts for the overalp of the WL states with the Mn, y = (3^ is the chemical potential 
of the WL carriers, z = {Em — E'^i)/kBTh is the change in energy of the Mn spin divided 



for 



by the effective temperature of the WL carriers, and I{y, z) = kBTh-^^y—^Log 
z > Q and I{y, z) = el^l/(?/, —z) for 2; < 0. 

The solution of the master equation, taking as initial condition p+5/2 = 1, is shown in 
figure (I7j). In the left panel we plot the evolution of the occupations of the 6 Mn spin 
projections along the growth axis. The decay of p+5/2 is accompained by a rise of the 
+3/2 state. This state is also connected with the 1/2 state whose occupation starts to 
build. Interestingly, the population of the +3/2 states overshoots initially. The figure 
shows a "falling domino" effect. They evolve towards thermal equilibrium. The average 
magnetization {M)(t) = J2n=i 6 ^^iPnit) decays according to 

(M)(t) = Moe-*/^"^ (6) 
16 



in all the simulations performed. In the figure we take Th — lO-ft', e = 0.1,?7 = l,P = 7jj,eV, 
m*=0.4, L = 20A, and rih = lO^^cm^^ and we obtain Tim — 8 ns. Thus, this mechanism 
is consistent with the Mn spin relaxation time that we observe. In the figure ([7]) we see 
how Tim increases as the density of carriers decreases. The order of magnitude of Tim is 
consistent with the calculations for Mn spin relaxation induced by scattering with electrons 
in quantum wells^. There Tim — 2 x 10^ ns is obtained in a quantum well with L = 80A. 
Taking into account that the electron Mn coupling a is 4 times smaller than (3 and the L is 
4 times bigger, there Tim for holes should be at least 16^ shorter. 




FIG. 7: (Color online), (a) Pn for Mn spin states. These can be considered as the conditional 
probabilities that, having observed the Mn spin at t = in the state = +5/2, a state with 
spin Mz is observed at time t. (b) Tim as a function of hole density (see text), (c) Average 
magnetization corresponding to panel (a), in logarithmic scale. 

Thus Korringa exchange with a sufficiently large density of photoholes might account for 
the Mn spin relaxation when the dot is in the optical ground state. We consider now a 
different mechanism, involving single hole charging of the dot, motivated by the fact that 
the PL often shows a weak contribution of positive trions. 

2. Mn spin relaxation due to charge fluctuation in the dot 

The elementary process for spin relaxation due to charge fluctuation is a combination of 
3 steps. First, a hole tunnels into the empty dot with the Mn in the spin state M. Second, 
the hole and the Mn exchange spin, so that the Mn is now in the state Mil. Third, the 
hole tunnels out of the dot. Thus, in the master equation we need to include both the 6 
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charge neutral states of the Mn doped dot as well as the 12 states with one hole inside. 
The dissipative dynamics, induced by V connects neutral states with charged states (+) and 
viceversa (— ). The charging transition rates are given by 

rW = ToriF {{E^ - Err.))\{n\ |0)(^||m)p (7) 
and the decharging rate: 

= To (1 - np{E^ - EJ) \{n\ \0)Mm)\' (8) 

cr 

Here npix) is the Fermi function which depends on the temperature and chemical potential 
of the electron reservoir, En and Em are the eigenstates of the QD Hamiltonian and Tq = 
Ylk I^Pp(-^/) tunneling rate for the fermion in and out of the dot. The matrix 

elements featured in eq. ( ITIIHI) are strongly spin dependent and only connect states in which 
the Mn spin changes, at most, by one unit. 

An important quantity in our simulations is the average charge, (g) = J2n=7 laPni'^)- 
Since we consider an initally neutral dot, this quantity goes from zero to the steady state 
occupation in a time scale given by Tig = (Fq)"^. In this second set of simulations we 
measure time in units of Tig. The evolution of the magnetization and the empty dot occu- 
pations on this case look very similar to those of figure (^^. In particular, the magnetization 
decays exponentially in a time scale much longer than the charge relaxation time. 

We have computed the Mn spin relaxation time, Tim, as a function of the parameters 
of the simulation, Fq, fc^T, J, e and (q). We have found that the crucial quantities are the 
average charge, the spin flip anisotropy th and Fq, as shown in figure ([8]). To very good 
approximation, we have that FqTim depends linearly on {q)~^, for moderate values of (q). 
Thus, this mechanism shows a strong dependence on the density of carriers in the WL, in 
agreement with our observations. 

In figure [8] we see that, in the case of tunneling electrons (e = 1), for an average charge 
of 0.1 we have Tim — IOOFq^. For tunneling holes with eh = 0.2 and (g) = 0.1 we have 
Tim — 600Fq ^. Interestingly, these numbers are weakly dependent of the value of J and on 
the temperature. In order to obtain an absolute number for T^^ we need an estimate of Fq 
the charge scattering rate. An upper limit for Fq is provided by the linewidth observed in 
the PL spectrum. Since the charge scattering would result in a broadening of the linewidth, 
we can infer that HTq < 50/ieV. Thus, using h = 0.65meVxps we have T*^ > 600^ ~ 8ns 
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FIG. 8: (Color online). Mn spin relaxation time, originated by charge fluctuations in the dot, as a 
function of average charge in the dot for electrons and holes with different eh values. The straight 
lines are linear fits of the numerical data. 

for (g) = 0.1. This is a lower bound for the spin scattering time, or an upper bound for the 
spin scattering rate. Thus, single hole tunneling events could relax the spin of the Mn in 
a time scale of 8ns for a dot which is charged 10 percent of the time. Agreement with the 
experimental result can be obtained by taking smaller average charge in the dot or smaller 
spin-flip interaction or smaller tunneling rate Fq. We note that the proposed mechanism is 
similar to the hole spin-relaxation mechanism proposed by Smith et alM. There, a spin up 
hole tunnels out of the dot and spin down tunnels in, resulting in an effective spin relaxation 
for the hole. 

In summary, we have two mechanisms that account for carrier induced Mn relaxation in 
the dot in a time scale of 10 nanoseconds. In both cases the dot must be considered as an 
open system. In the Korringa relaxation mechanism the Mn spin is exposed to extended WL 
states. In the charge-fluctuation mechanism the occupation of the dot state fluctuates, due 
to its coupling to extended WL states. These two mechanisms are not exclusive and might 
operate at the same time. Whereas the Korringa mechanism is almost identical to that of 
(Cd,Mn)Te quantum wells^«^ the charge- fluctuation mechanism is specific of quantum dots 
coupled to a continuum. 
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VIII. CONCLUSION. 



In conclusion, we used time resolved and photon correlation measurements to probe the 
spin dynamics of a single magnetic atom (Mn) interacting with photo-created carriers at zero 
magnetic field. Fluctuations of the localized Mn spin control the statistics of the photons 
emitted by a single Mn doped QD. The spin relaxation time of the Mn atom, around 20ns, 
appears to be particularly sensitive to the spin-spin coupling with free carriers injected 
in the QD or on the vicinity of the QD. The fast photo-induced spin relaxation of the 
Mn is dominated by spin and energy transfer from hot photo-carriers directly to the Mn 
atom. Minimizing the injection of free carriers by optimizing the excitation conditions would 
make it possible to reach longer fiuctuation times. As a consequence, a determination of the 
intrinsic Mn spin dynamics at zero magnetic field, ultimately limited by the spin-lattice 
interaction, could only be performed in the absence of free or confined carriers. 

This work was supported by French ANR contracts MOMES and CoSin. 



Electronic address: lucien . b esomb es@grenoble.cnrs.fr 



^ M. Xiao, I. Martin,E. Yablonovitch,H. W. Jiang, Nature (London) 430, 435 (2004). 

2 J. Fernandez-Rossier and R. Aguado, Phys. Rev. Lett. 98, 106805 (2007) 

3 J. Lambe and C. Kikuchi, Physical Review, 119, 1256 (1960). 

^ M.K. Kneip, D. R. Yakovlev, M. Bayer, A. A. Maksimov, I. I. Tartakovskii, D. Keller, W. 
Ossau, L. W. Molenkamp, A. Waag, Phys. Rev. B 73, 045305 (2006) 

5 T. Dietl, P. Peyla, W. Grieshaber, Y. Merle d'Aubigne, Phys. Rev. Lett. 74, 474 (1995). 

6 A.J. Heinrich, J.A. Gupta, CP. Lutz, D.M. Eigler, Science 306, 466 (2004) 

^ D. Kitchen, A. Richardella, J.-M. Tang, M. E. Flatte, and A. Yazdani, Nature (London) 442, 
436 (2006). 

^ L. Besombes, Y. Leger, L. Maingault, D. Ferrand, H. Mariette, J. Cibert, Phys. Rev. Lett. 93, 
207403 (2004) 

9 L. Maingault, L. Besombes, Y. Leger, C. Bougerol, H. Mariette , Appl. Phys. Lett. 89, 193109 
(2006). 

^° C. Couteau, S. Moehl, F. Tinjod, J.M. Gerard, K. Kheng, H. Mariette, J. A. Gaj, R. Romestain, 



20 



J.P. Poizat, Appl. Phys. Lett. 85, 6251 (2004). 

Y. Leger, L. Besombes, L. Maingault, D. Ferrand, H. Mariette, Phys. Rev. B 72, 241309 (2005). 
L. Besombes ,Y. Leger, L. Maingault, D. Ferrand, H. Mariette, J. Cibert, Phys. Rev. B 73, 
161307 (2005). 

^2 J.M. Smith P. A. Dalgarno, R.J. Warburton, A.O. Govorov, K. Karrai, B. D. Gerardot, P. M. 

Petroff, Phys. Rev. Lett. 94, 197402 (2005) 
14 A. Kiraz, S. Flth, C. Becher, B. Gayral, W.V. Schoenfeld, P. M. Petroff, Lidong Zhang, E. Hu, 

A. Imamoglu, Phys. Rev. B 65, 161303(R) (2002) 
1^ D. Scalbert, J. Cernogora, C. Benoit la Guillaume, Solid State Commun.66, 571 (1988). 
16 A. O. Govorov, A. V. Kalameitsev, Phys. Rev. B 71, 035338 (2005). 
1^ D.H. Feng, LA. Akimov, F.Henneberger, Phys. Rev. Lett. 99, 036604 (2007). 
1^ A. Vasanelh, R. Ferreira, and G. Bastard, Phys. Rev. Lett. 89, 216804 (2002) 
19 M.G. Tyazhlov, A.I. FilHn, A.V. Larionov, V.D. Kulakovskii, D.R. Yakovlev, A. Waag, G. 

Landwehr, JETP 85, 784 (1997). 
2° L. Besombes, Y.Leger, L. Maingault, D.Ferrand, H.Mariette, Phys. Rev. B 71, 161307(R) 

(2005). 

21 L.M. Woods, T. L. Reinecke, and R. Kotlyar, Phys. Rev. B 69, 125330 (2004). 

22 Y. Lcgcr, L. Besombes, L. Maingault, H. Mariette, Phys. Rev. B 76, 045331 (2007). 

23 J. K. Furdyna, J. Appl. Phys 64 R29 (1988). 

24 J. Lambe and C. Kikuchi, Phys. Rev. 119, 1256 (1960). D. Scalbert et al, Solid State Commu- 
nications 66, 571 (1988). 

2^ B. Koenig, I. A. Merkulov, D. R. Yakovlev, W. Ossau, S. M. Ryabchenko, M. Kutrowski, T. 

Wojtowicz, G. Karczewski, and J. Kossut Phys. Rev. B 61, 16870 (2000) 
26 E. Souto E, O. A. C. Nunes , A. L. A. Fonseca, Solid State Comm. 129, 605 (2004) 
2^ Y. Leger L. Besombes, J. Fernandez-Rossier, L. Maingault, H. Mariette , Phys. Rev. Lett. 97, 

(2006) 

28 J. Fernandez-Rossier, Phys. Rev. B.73, 045301 (2006) 

29 J. Fernandez-Rossier and L. Brey, Phys. Rev. Lett. 93 117201 (2004) 

3° Claude Cohen- Tannoudji, Jacques Dupont-Roc, Gilbert Grynberg Atom-Photon Interactions: 

Basic Processes and Applications, John Wiley and Sons, New York, 1992. 
31 M. Qazzaz, G. Yang, S. H. Xin, L. Montes, H. Luo, J. K. Furdyna Solid State Comm. 96, 405 

21 



(1995) 

J. M. Smith, P. A. Dalgarno, R. J. Warburton, A. O. Govorov, K. Karrai, B. D. Gerardot, and 
P. M. PetrofF Phys. Rev. Lett. 94, 197402 (2005) 

P. A. Dalgarno, M. Ediger, B. D. Gerardot, J. M. Smith, S. Seidl, M. Kroner, K. Karrai, P. M. 
Petroff, A. O. Govorov, and R. J. Warburton Phys. Rev. Lett. 100, 176801 (2008) 



22 



